Eur ophysics Letters PREPRINT 



Boltzmann conjecture, meta-equilibrium entropy, second law, 
chaos and irreversibility for many body systems. 

PlERO ClPRIANl(*) 

MIUR - (Italia) 

ICRA - P.le Repubblica, 10-12, 65122 Pescara (Italia) 



PACS. 05.20.-y - Classical statistical mechanics. 

PACS. 05.70.Ln - Nonequilibrium and irreversible thermodynamics. 

PACS. 45.05.+X - General theory of classical mechanics of discrete systems. 



Abstract. - A heuristic generalization of the Boltzmann- Gibbs microcanonical entropy is 
proposed, able to describe meta-equilibrium features and evolution of macroscopic systems. 
Despite its simple-minded derivation, such a function of collective parameters characterizing 
the microscopic state of N-body systems, yields, at one time, a statistical interpretation of 
dynamic evolution, and dynamic insights on the basic assumption of statistical mechanics. Its 
natural (implicit) time dependence, entails a Second Law-like behaviour and allows moreover, 
to perform an elementary test of the Loschmidt reversibility objection, pointing out the crucial 
relevance of Chaos in setting up effective (statistico-mechanical and dynamical) arrows of time. 
Several concrete (analytical and numerical) applications illustrate its properties. 



Introduction. - The definition of an entropy-like quantity in terms of dynamical vari- 
ables for mechanical systems out of equilibrium, is a fundamental issue of Statistical Mechanics 
(SM). In the last decade, most studies focused on Non Equilibrium Stationary States (NESS) 
of driven systems, [1]. In parallel, however, there has been a revival, [2] of the investigations 
on the seemingly simpler issue, dealing with isolated systems, left free to evolve from a non 
equilibrium initial state; i.e., on the original concern of Boltzmann: the dynamical interpre- 
tation of the tendence of macroscopic systems to spontaneously relax to an equilibrium state, 
and why and how natural evolution is associated with an increase of a macroscopic state vari- 
able, called Entropy in the realm of Thermodynamics (TD), for which a definition in terms of 
microscopic phase variables is sought [3]. In these latter papers in particular, a modification 
of the (perfect gas) H function is suggested, when the contribution of the interaction poten- 
tial to the total energy is not negligible. These works acknowledge the deep (and generally 
ignored) Jaynes critical analysis, [4], on the (ab)use of the identification of the Boltzmann-H 
function with the (negative of the) thermodynamic entropy, and on the (logically equivalent) 
statement that the derivation of any flavour of a H-theorem should amount to a proof of the 
Second Law. 

Here, I present a further attempt, inspired again by the Jaynes thought-provoking analysis, 
which, far from giving a mathematically rigorous derivation, relies upon a naive extension 
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of a well known (and implicitly commonly adopted) generalization of the Boltzmann-Gibbs 
(BG) microcanonical entropy and does not require any modification of the most simple and 
fundamental mechanical definition of entropy, [5]. The interconnections with other recent 
works around this area, either from the viewpoint of an axiomatic treatment of t he Second 



Law, [6], or dealing with NESS and related chaotic hypotheses, [1,7] will appea af 1 )! elsewhere, 
[8], where conceptual, analytical and numerical aspects are also discussed in more details. 

In the following, I describe the basically simple idea, and how it, leaving formally un- 
changed the BG entropy, naturally leads to an implicitly time-dependent quantity. Then I 
apply it to three different models, discussing a) the analytic predictions; b) the surprising 
agreement between these and the numerical outcomes, in equilibrium and non-equilibirum 
cases; c) the experimental increase of generalized entropy along dynamical trajectories and d) 
the evidence of the crucial relevance of Chaos for the onset of irreversibility. 

The basic assumption. - The minimally biased choice for a reconciliation between Dy- 
namics, spontaneous evolution and statistical description of isolated systems, has to be micro- 
canonical in spirit: I will deal with N-body systems (in a /—dimensional configuration space), 
governed by a Hamiltonian: 

N 2 N ~ 1 N 

H (£'E)=E|^ + E E V ij (\\c li - qj \\)=E = )C + U . (1) 

i—l i—1 j=i-\-l 

The SM over the (2fN — l)-dimensional constant energy surface, establishes the logical 
connection between Dynamics and TD, exploited through the BG entropy, Sbg(E, N,V) = 
ks In W(E, N, V), where W is the measure of T, E (k B being the Boltzmann's constant), and 
depends on few macroscopic parameters (E,N and, e.g., the spatial volume V). 

However, it is well known that, in the case of isolated systems, dynamical trajectories 
actually explore only a (zero-measure) subset of liyng down rather on the D = (2/iV — 1 — 
M)-dimensional manifold, fl, defined by E and the other M integrals of motion of the system. 
Indeed, if a consistent SM is sought, Sbg has to be suitably modified, whenever the dynamics 
admits additional conserved or costrained quantities. This generalization is simple and is 
tacitly assumed in most textbooks and discussions, mainly because, generally, it is M <C N and 
the results do not depend on the choice (see, however [10] for a curious counterexample). The 
above definitions lead to the microcanonical SM and work egregiously well for equilibrium TD. 
However, it has been argued (see, e.g., [11]) that the above expressions, if suitably interpreted, 
could give a measure of entropy even in noncquilibrium states. Although such a terminology 
will horrify anyone eager for rigorous definitions, Jaynes instead was a convinced supporter of 
this interpretation, though referring to the Gibbs entropy, Sq — — fcs J^phipdT, [4]. 

In the same spirit, the basic extension here proposed consists in assuming that the region f2 
whose measure, W = , enters in the definition of the generalized entropy, depends not only 
on the (conserved) values of energy, H(q,p) = E and on the additional integrals of motion, 
JZ(q,p) = I, where I = {I±, . . . ,Im}, fixed a priori by the initial conditions and thereafter 
rigorously conserved, but can have a functional dependence also on other collective variables, 
{Vi}, whose evolution is instead determined (apart from statistical fluctuations), implicitly 
and self-consistently, by the dynamics. As the {Vi} evolve, according to Hamilton's equations, 
the measure of the effectively available phase-space volume changes, so that the generalized 

( 1 ) This paper start to summarize studies and reflections, formed during my Ph.D. studies, [9]. The revivaJ of 
the investigations on the dynamical foundations of the Second Law and the wish to acknowledge the conceptual 
relevance of the work of E.T. Jaynes, persuaded me, eventually, to collect those old notes. 
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entropy function, Sbg, of the system is allowed to evolve as well: were 57 = fl(E, N, I) alone, 
then, clearly, no entropy evolution could be explained. 

The basic consistency requirement amounts to verify whether the quasi deterministic evo- 
lution of the {Vi} implies almost always an increase of Sbg, and, viceversa, whether the 
maximization of the latter with respect to the macroscopic parameters, gives predictions about 
their behaviour consistent with the dynamics. In the affirmative case, a meaningful (yet heuris- 
tic) reconciliation between Dynamics and TD is obtained. It has to be emphasized how, in 
this scheme, physical intuition plays a major role to single out the best suited macroscopic 
parameters able to describe the dynamical and statistical evolution of the system. All selec- 
tion criteria should consider that TD and SM, deal with macroscopic systems; consequently, 
a possibly singular behaviour of single microscopic constituents, which does not alter the val- 
ues of collective parameters, cannot modify the pertinent level of description. To show the 
effectiveness of the approach, I will confine to few (yet physically relevant) cases, where the 
above selec tion is more or less straightforward and the physical meaning of the parameters 



transparen t 1 Therefore, introducing the specific energy, e = E/N , and assuming that the 



{Vi} vary within suitable intervals ILi(s, N; I), the generalization suggested reads: 

S BG {e,N,l;{Vi}) = k B hxW = k B ]n f 5[E - H{q,p)]8( M \l-K{q,p)]dqdp. (2) 

J Vie rii 

The (implicit) time dependence is hidden in the dynamically determined evolution of the {Vi}. 
Such macroscopic variables can be often singled out analyzing the scaling properties of the 
model. A common example, exploited below, is the virial ratio, Q = Ay(K,)/(U), between 
the (time averages of) kinetic and potential energies. A purely dynamical theorem (see, e.g. 
[12, 13]) assures that it attains, for confined motions, a well defined value, dependent, in the 
general case, only on the form of potential V and, possibly, on e. The phase-space volume can 
often be expressed as a function of the related instantaneous virial ratio, Q(t) = A\;lC(t)/U(t). 

Applications. As a first example let us consider the simplest N-body hamiltonian 
system (except, perhaps, the perfect gas): a chain (/ = 1) of N harmonic oscillators with 
nearest neighbour interactions. Although this system is integrable, for irrational frequency 
ratios, kinematic phase mixing takes place and motions of single oscillators decorrelate. It 
is obvious that no irreversible evolution can occur; however, while some variables approach 
a definite value only in time averaged sense, others assume, even locally, values practically 
coincident with asymptotic ones. For example, in this case, Qho = (IC)/{L4) = 1 is energy 
independent and Q(t) oscillates around this value, within statistical 0(N~ 1 ^ 2 ) fluctuations, 
unless collective oscillations occur, which are, however, quickly damped out by the (kinematic) 
mixing. As a first exercise, [9], I show how the dynamic evidence, Q = 1, receives within the 
present approach a statistical explanation (or, at least, interpretation). Indeed, assuming 
for simplicity m, = m and introducing suitable scales, (R, V), for positions and momenta 
coordinates, q.; = i?Xj, and = mVyi, such that 2e = mV 2 + mu> 2 R 2 , with u> 2 a weighted 
average frequency and exploiting the definition of Q(t), it is possible to write, [8], 

v>= ? Q ",;R>= - 2£ , (3) 

m(l+Q) mtu 2 (l + Q) 



( 2 )In particular, dealing with self-confining interparticle potentials allows to neglect the effects of a vessel 
and to exploit the fact that, once E and N are fixed, the spatial volume V is implicitly costrained within well 
defined limits, at least on timescales of interest, see [8]. 
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from which, substituting into eqs.©, we obtain 

4Qe 2 



5 f N L 1 



w 2 (l + Q) : 



■ y) 



(4) 



with I(x, y) a dimensionless integral, depending only on the microscopic state, and evolv- 
ing (much) more slowly than, Q. Despite its simplicity, eq.(|3]) has two important conse- 
quences: a) the maximization of Sbg with respect to <2, yields the virial equilibrium condition: 
OSso/dQ = ^==> Q = 1; b) the microcanonical temperature, T = (dSsG/dE) -1 , leads 
to the well known formula for the heat capacity of a (1-D) harmonic crystal: E = NksT, 
coherent, incidentally, with the orthodicity requirement, [15]. 



Before to go beyond what could have been considered just a funny educational curiosity^ 
I want to point out that, in the general case, the evolution of the collective parameters is not 
the only source of entropy increase; however, a separation of scales (of size and time) exists, 
such that the gross and faster evolution is driven by these quantities, [8, 14]. 

Let us consider now a system (called here RMLB-model) which has been used to model 
phase transitions in small clusters of atoms, [16], whose Hamiltonian is 

N 2 N 



H(q,p) = ^L. + ^K Q Nr 2 [l - br 2 + cr 4 ] , with r 2 = A^ 1 ^ q 2 . (5) 

i— 1 1 i—1 



The TD (and, surprisingly, the dynamics too) of this model presents (even in the case Kq = 1, 
c = 1 and m. L = 1) several interesting aspects, that have been, [16], or will be discussed 
elsewhere in details, [8]. The approach proposed, allows to recover its TD peculiarites and non 
trivial phenomenology (as, e.g., multiple local entropy maxima and the occurrence of negative 
heat capacities) and to interpret the connections between dynamical and TD computations, 
even in a system, again, without any dynamic stochasticity. 

Assuming, as before, / = 1, exploiting, in this case, the virial theorem, [8], and introducing 
the specific entropy, o = Ssc/NkB, we get a parametric entropy- energy relation through the 
order parameter p = (r 2 ) 1 / 2 : 

e = (P 2 - | P 4 + 2p 6 ) ; a S \ In [2p 2 (p 2 - 26/ + 3p 6 )] , (6) 
As before, the microcanonical Temperature, which reads here: ksT(e) = 



(with 

(da/dp) p=pe 

p x such that s(p x ) = x)), yields the caloric curves, and confirms also for this model, that, 
within this approach, the orthodic condition holds, [8]. To compare concretely the predictions 
of the proposed SM with the outcomes of Dynamics, I've performed several Molecular Dy- 
namics (MD) simulations of Hamiltonian §5§ for different values of b and e and computed the 
corresponding dynamical temperatures, through well known formulas, [17]. These are com- 
pared, in figlU with the theoretic curves derived above. We see that the agreement is absolute 
almost everywhere, except in the regions of negative heat capacity. As anticipated above, the 
interesting features of this model, indeed, originate from the existence of a threshold value, 
b > 6*, above which coexistence regions appear. There, as shown in the insets of figlU both 
e(p) and o~(p) are not anymore invertible functions, so that p e is multivalued as is T{e). This 
is not at all a drawback: on the contrary, the big fluctuations there occurring, help to un- 
derstand the nature of the coexistence in this system: in a (microcanonical) MD simulation 

( 3 ) And I've been so persuaded fifteen years ago, even if it was not the only example I had found at that time. 
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Fig. 2 - Comparison of the theoretical, microcanonical distribution of the Q-values, predicted from 
eq.(0, with the experimental, relaxed distribution, computed along a numerically integrated hamil- 
tonian trajectory, for a SGS with TV = 1000. The dashed curve represents a gaussian fit. In the inset 
the same comparison is made, for N = 16000 but for a much shorter time. 



e is fixed, though this doesn't fix obviously r. In a pure phase, b < &», however, it oscillates 
in a single range, so that p attains a well defined value; if b > 6*, instead, r can be trapped, 
intermittently, within different bands, compatible with the given e. Correspondingly, p and, 
consequently, cr and T, show much larger fluctuations. 

The dynamics of both models analyzed so far is not enough unpredictable, and this makes 
them not exactly the kind of system one expects to be suited for a SM approach. Despite 
this, the results shown allow to argue that deterministic chaos is not always essential for a (at 
least partially) successful SM description of (meta-)equilibrium states. 

To investigate deeply the distinguishing features of TD, e.g. the Second Law and the 
origin of an effectively irreversible evolution, I will focus now on the (3D) self-gravitating 
N-body system (SGS), whose dynamics is surely unpredictable, has been the original stim- 
ulus which brought about, [9], the first reflections and for which preliminary accounts have 
already appeared, [14]. For simplicity, I assume to have a boun d (e = E/N < 0) system 
of N equal mass (m, = m) particles, interacting via a softene d ( )l Newtonian potential, 



( 4 ) As discussed elsewhere, [8], the softening is adopted here not (only) to avoid numerical singularities, 
but mainly to be fully coherent with the basic assumptions of the approach proposed. Therefore, in what 
follows dij = *J (rj — r j)' 2 + n 2 , where 77 = r?* d. That is, the two-body interaction is softened when the 
particles get closer than a (small) fraction of the mean interparticle separation, d = GN 2 / 3 m 2 /\e\. In all the 
numerical simulations here discussed, every particle has unit mass, the initial conditions are generated at virial 
equilibrium, — 2/C = U, with specific energy fixed e = —10, and particles distributed initially with uniform 
spatial density in a spherical box of radius Rq, with a velocity distribution approximately maxwellian. Time 
is measured in units of dynamical (or crossing) time [9,14,18], to = (Gg)~ 1 / 2 = N/\s\ 3 ^ 2 . The numerical 
integrations are performed in double precision, using a high-order symplectic algorithm, [20] , and the time-step 
is chosen small enough to keep the (maximum) relative energy error below 10 -4 . The position and the velocity 
of the center of mass of the system are rigorously conserved within the round-off errors, i.e., within 10 -14 . 
The consequences (or lack thereof) of modifying some of the choices are discussed elsewhere, [8]. 
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Vy(||qj — qyH) = —Gmirrij j d%j . Gravitational interaction is not rigorously confining] it can be 
argued, however, [9,14,18] that, up to any physically meaningful time, the number of escaping 
particles is negligible. Thus, as before, we can introduce suitable scale factors for coordinates 
and momenta, defining the gravitational and inertial radii, Rg = QT, mi) 2 / Q^^y trumj / dij) ~ 

1 /2 

GNm 2 /\e\ and Ri = IT^ m^ 2 / ^2 m i] , respectively, and the ratio among them, a = Ri/Rq. 
After some elementary algebra, [8], we can express again the phase-space volume as a func- 
tion of standard parameters (e,N), collective phase functions, Q = —2IC/U and a, and a 
dimensionless integral, J , accounting for the correlations amongst particles, to obtain: 



~ _ 3Nk B 

o bg = Nkbct = — - — in 



G 2 N 2 m 5 - , ^ x . / 
a 2 Q(2-Q) J(x,y) 



(7) 



The different kinds of parameters in eq.([7]) play distinct roles: as far as we consider the system 
isolated, e and N (and clearly m) are rigorously constant, whereas Q, a and J evolve, though 
on very well separated (and increasing) timescales. A thoroughful analysis of their hierarchy 
is presented elsewhere, [8,14]. The implications of eq.® follow immediately: a) the definition 
of the microcanonical temperature leads to the caloric curve for SGS, 

T- 1 = {dS BG /dE) = k B (da/de) => e = -h B T , (8) 

yielding a direct proof of the well known (for SGS) negative heat capacity, b) As in pre- 
vious cases, Sbg m eq.© is maximized for Q = 1, and this confirms again the statistical 
intepretation of virial theorem, c) What's more, for this system, the present approach allows 
to go well beyond the mere prediction of the correct average value: eq.® provides also a 
distribution law for the values of Q in the microcanonical ensemble. This is a rather con- 
vincing result if we ponder over the surely chaotic nature of SGS, [9,14,19]; implying that 
their statistical properties are likely to be very strong. It is therefore reasonable to expect 
that (microcanonical) ensemble probability distributions should be accurately sampled during 
a long enough hamiltonian trajectory. The agreement between Dynamics and (modified) mi- 
crocanonical SM is clearly evident from Figure [2 the histogram reports the effective counting 
of the actual values assumed by Q(t) along a numerically integrated trajectory (N = 1000 
and 77* = 0.01)) while the solid line represents (up to a normalization factor) the pdf for 
the Q- values, as predicted by eq.©. It can be argued that both distributions are very close 
to the gaussian curve (dashed line), which seems to reproduce equally well the numerical 
data. However, there is a physically important distinction: the gaussian results from an opti- 
mized empirical fit, whereas the /^-canonical pdf has no fitting parameters. The difference is 
strinkigly evident in the inset, where an analogous comparison is made for a simulation with 
N = 16000, but on a much shorter time, just comparable with the virialization time. There, 
the gaussian fit is apparently unsatisfactory, as the weight of the (relatively) large oscillations 
of Q during the virialization process, imposes to the fit a very large width. On the other 
hand, it should be stressed that, on increasing N (and awaiting for a long enough integration 
time), the ^(-canonical pdf and the gaussian fit will overlap, as the former, in the N —> 00, 
approach a gaussian with width a 2 oc iV^ 1 and the amplitude of the fluctuations will vanish, 
so that the gaussian fit shrinks, [8]. Thus, the proposed framework, allows to reproduce the 
statistical distribution of collective parameters along hamiltonian trajectories even for finite 
time, when equilibrium distributions fail, recovering, however, in the TV — > 00 and t/tfj —> 00 
limits, the normal equilibrium results. So far we have thus seen that a maximum entropy-like 
criterion applied to Sbg is n °t only consistent with dynamical equilibrium conditions but 
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Fig. 3 Fig. 4 

Fig. 3 - Increase with time (measured in units of the dynamical time, tr>) of a* , along hamiltonian 
trajectories of several SGS. From top to bottom it is N = 1000, 2000, 4000, 8000, 16000. Except for 
the largest N, the plots represent the average of different realizations. Always it is n* = 0.02. 

Fig. 4 - Reversibility objection and Chaos. Evolution of <r* along single trajectories of N = 4000 
SGS's with n* = 0.02, 0.2, 2, from top to bottom. At t = Unv = 20to the velocities are reversed. The 
upper left inset is an enlargement around ti nv for the smallest n*; the other inset contains the plots 
of 8*(t) = a* {t) - a*(2t inv - t). 



predicts also correctly (meta-)equilibrium statistical and TD properties of many body sys- 
tems. However noteworthy may be that, one could wonder whether Sbg satisfies also deeper 
properties expected to hold for an entropy function. If we were able to show (if not to prove!) 
that Sbg increases almost monotonously with time, then we could be really confident that 
the naive estimate of the phase space volume £l(E, N, I; {Vi}) contains a profound meaning, 
connected with the Second Law and the issue of irreversibility. That this is indeed the case 
is, I thi nk, unambiguously shown in figure [31 where the (numerical) time evolution of the 



quantit y 5 )l a* = In [a 2 Q(2 — Q)l, for different values of N. After a relatively fast (with a 



iV-independent rate) initial increase, a* steadily grows, although with a slower rate the larger 



is TV, [8], and with fluctuations which also are smaller the greater is M( k Sbg, computed 
along (microcanonical) dynamical trajectories behaves like the TD entropy and, moreover, 
the frequency of local (in time) violations of the Second Law vanishes for large N. 

This is much of though not all, the content of the Second Law: this law establishes indeed 
what is (popularly) known as an arrow of Time; that is, it enables to distinguish the past 
from the future evolution. This contrasts with Hamiltonian dynamics, invariant under time 
reversal. I will not discuss here [8] the many attempts to settle the issue, but simply recall 
two amongst the many arguments invoked to justify the onset of macroscopic irreversibility in 
the dynamics of mechanical systems: the occurrence of non typical initial conditions and the 
presence of a chaotic dynamics. Indeed, discussing the first two models, we have seen that 

( 5 ) Which is, essentially, the only evolving contribution to the specific entropy Sbg/^^b, as the other time- 
dependent quantity, the integral J , evolves on much longer timescales. Notice that time is measured in units 
of the so-called dynamical time, to oc TV. 

( 6 ) The details and the implications in astrophysical context of these results will be discussed elsewhere, along 
with the clear-cut intensive nature of o~* emerging from the figure. 
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even systems which are definitely not chaotic, possess (some) collective properties which can 
be described statistically. This amounts to say, that TD concepts can be useful to describe at 
least some equilibrium properties of even dynamically regular macroscopic systems, provided 
that no reference is made to non equilibrium behaviour, associated to the Second Law. 

With regard to the issue of irreversibility, were the entropy increase due solely to the 
peculiarity of the initial state, then the Loschmidt reversibility objection would not lead to 
any paradox, and entropy could possibly decrease going back in time. I claim, however, 
that this is not the case, and precisely that Chaos plays an essential role in the true onset 
of irreversibility. In this perspective, therefore, the relationship between Second Law and 
irreversibility loses a tiny part of its complete and absolute logical equivalence: the first can 
originate from a non-generic initial state of the system, the impossibility to retrace exactly 
the evolutive path backward in time arises from the presence of a chaotic dynamics. In this 
evenience, consequently, entropy functions exist that increases in both directions of time, 
irrespectively of the choice of initial conditions. 

Figure 2] unequivocally supports the above assertions, giving them a clear empirical ev- 
idence. The curves there shown represent, as in figj3j the time evolution of a* for three 
N = 4000 particles systems, all starting from the same initial conditions, but with vary- 
ing softening parameter, n. A t va riance with the previous figure, however, the validity of 
Loschmidt reversibility obiectio A( 7 )\ is checked, reversing, at time t — t inv , every particle's ve- 
locity. The results are clear: A) for a strongly, or even moderately, chaotic system, the entropy 
evolution is reversed only for very short time intervals after the inversion; subsequently the 
(dynamical) instability drives again the collective variable a* to increase, obeying to a Second 
Law even for the backward time evolution (see upper inset). B) Viceversa, if the dynamics 
is virtually regular, as in the case of a very large 77*, although an initial entropy increase oc- 
curs as well, due to the non complete genericity of the initial state, after the inversion, the 
macroscopic quantity a*, traces back along its path, taking exactly the same values assumed 
prior to inversion, and accompanies the return of the system to the initial (macroscopic and 
microscopic) state (lower inset). 

I'm obviously aware that the above results deserve a deeper discussion [8] and that the 
interpretation of the numerical effects mimicking the unavoidable external disturbances should 
be carefully investigated; nevertheless, I find it extremely instructive, as it shows a clear link 
between microscopic instability (let it be driven by numerical or ambiental noise) and irre- 
versible behaviour of macroscopic quantities, whether it is agreed upon to call them entropies 
or not. 

Conclusion. - A simple naive estimate of phase space volume, compatible with the values 
of macroscopic parameters of several hamiltonian N-body systems has been used to introduce 
a simple generalization of the BG entropy that obeys and satisfies the corresponding Second 
Law. The results here presented are unambiguous and do not require anything more than an 
undergraduate mathematics. The same conceptual approach can (and will [8]) be generalized 
to a wider class of systems, just using a more refined formal treatment. Admittedly, the 
extreme simplicity and the heuristic nature of the approach call perhaps for a more rigorous 
formalization. Nevertheless, its relevance resides just in its elementary derivation and dis- 
closes itself through the ability to describe adaptively the TD properties of a system in terms 
of suitable macroscopic parameters which evolve on suitable timescales. Another conceptually 

( 7 ) It is wise to recall that Loschmidt's objections (as well as Zermelo's one) were against the Boltzmann's H- 
theorem, so related to the interpretation of the Boltzmann's _ff-function, based on the single particle distribution 
function. Nevertheless, it is clear that Loschmidt's criticisms could apply equally to the present context 
(whereas for the Zermelo's objection it is not exactly equivalent, [8]). 
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appealing feature of the approach here proposed is that the hypotheses needed for different 
levels of statistical descriptions emerged naturally, along with the development of the reason- 
ing. In particular, we have seen how the conditions, like literal ergodicity or mixing, usually 
required, are not always essential to justify an effective (though partial) TD description. 
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